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We generalize Dykhne's calculation of the effective resistance of a 2D two-component medium 
to the case of frictional drag between the two parallel two-component layers. The resulting exact 
expression for the effective transresistance, p^g-, is analyzed in the limits when the resistances and 
transresistances of the constituting components are strongly different - situation generic for the 
vicinity of the classical (percolative) metal-insulator transition (MIT). On the basis of this analysis 
we conclude that the evolution of p q across the MIT is determined by the type of correlation 
between the components, constituting the 2D layers. Depending on this correlation, in the case of 
two electron layers, p°g- changes either monotonically or exhibits a sharp maximum. For electron- 
hole layers p°g- is negative and |p°£j| exhibits a sharp minimum at the MIT. 



Introduction. Frictional drag between two layers has been first predicted theoretically [1,2] and later observed exper- 
imentally [3,4]. The characteristics measured in experiment is the drag resistance, p D = V p /I a , where, V p , is the 
voltage, built up in the passive layer upon passing the current, I a , through the active layer. Experimental observa- 
tions [3,4] have inspired a great number of theoretical studies of the frictional drag for different realizations of the 2D 
electron (hole) systems, constituting active and passive layers [5-20]. In parallel, a general formalism for calculating 
of drag was advanced [21-27]. 

In all theoretical papers published by now, the parallel layers were assumed perfectly homogeneous on the macro- 
scopic scales (usually, the scales exceeding the carrier mean free path, I). Incorporating tunneling bridges [13] or 
assuming correlations between the wave functions of the two layers [14] did not violate their macroscopic homogene- 
ity. Also, except for Refs. [17,18], the temperature was assumed high enough, thus allowing to neglect the mesoscopic 
fluctuations due to coherence of different regions of the layers. The question about the magnitude of drag between the 
2D layers, which are strongly inhomogeneous macroscopically was not addressed in the theories [5-27]. This question 
is studied in the present paper. 

To be specific, consider first the following situation. Assume that the passive layer is a good metal, kpl p 3> 1, which 
is perfectly homogeneous with a fixed concentration of electrons, n p = kp/2ir. The concentration, n a (r), of electrons 
in the active layer is determined by the concentration of donors, Nr>(r). Due to, say, imperfections in the doping 
process, Noir) fluctuates on a macroscopic scale with very large correlation length, R c ^> I (such an assumption was 
previously adopted in Refs. [28,29]). Assume now, that the active layer can be depleted by applying the gate voltage, 
V g . Without the gate voltage, V g — 0, we have n a (r) = Nn(r). Upon increasing V g , the electron concentration 
changes as n a {r) = No{r) — aV g , where the dimensionless coefficient a describes the depletion rate. Assume also, 
that at V g = the concentration, n a , is high enough, so that even with fluctuations, every region of the active layer is 
metallic. As V g increases, the local resistivities p a (r) = p a {Nrj(r) — ctV g } will also increase, but at different rate, so 
that within a certain domain of V g the inhomogeneities in Nr>(r) will become important. Namely, while some regions 
of the active layer will remain metallic with p a = p\ weakly dependent on temperature, the remaining area of the 
active layer will turn into insulator with p a — p^ oc exp(W/T), where U is the activation energy. Then it is clear, 
that at certain critical V g — Vg, the metallic regions will occupy exactly 50% of the area of the active layer. In other 
words, the classical MIT will take place within the narrow interval \V g — V g \ -C V g . The width, 5V g , of this interval 
can be related to the critical exponent, t ~ 1.3, of conductivity in the classical percolation [30]. Indeed, in the limit 
P2 — ► oc, the resistivity near V g diverges as p a {V g ) ~ pi [(V g — V g j /V g ] . Conversely, in the limit pi — > 0, we have 

P a (Vg) ~ Pi [(V g - V g c ) /V g c ]\ Then SV g is determined by matching the two behaviors, i.e. SV g /V g c = (pi/ ' p%) 1/2i \ 
For activated character of transport in insulator, assumed above, SV g shrinks with temperature as exp(— UjltT). 

It is important to note that, in addition to the above qualitative picture, there exists a sound quantitative result 
concerning resistivity at 2D classical MIT. Namely, as it was demonstrated by Dykhne [31], the exact value of p a at 

V g = Vg is equal to p a (Vg) = (plp?,) 1 ^ 2 - Moreover, the product p a (Vg + v) p a (V g - v) is equal to (plp?,) 1 ^ 2 for any 
v. Then the question arises about the behavior of the drag resistance, p^, in the vicinity of the classical MIT. 

It is obvious that, outside the interval \V g — Vg \ < Vg (pl/p^) 1 ^ 21 the effective transresistance is equal to p° on the 
"metallic" side of MIT and to p% on the "insulating" side, where p°, p% are the transresistivities between the regions 
with resistance pi and pf of the active layer and the metal of the passive layer, respectively. This is because, outside 
the transition region, the transport is dominated by the current paths going exclusively through the regions of either 
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low (metallic side) or high (insulating side) resistance. The main message of the present paper is that, similar to the 
value of p a (Vg) , the exact value of p D (Vg) = p°^ can be found. In particular, in the limit p°g. <C p^g, this value is 
given by 



To analyze the temperature dependence of p°^ one can use for pf a conventional expression for drag between two 
metals. Concerning the drag resistivity, p 2 , we have assumed that the transport in the insulating regions of the active 
layer is due to activated electrons. For these electrons, collisions with electrons in the passive layer, can be viewed as 
an additional source of scattering. From here we conclude, that both the conductance and trans conductance for the 
insulating regions are oc exp {—U/T). In transresistance, however, this exponent cancels out, so that the T-dependence 
of p 2 is weak. It is obvious from Eq. (1) that the magnitude of p^g- lies between pf and pf S> pf . Since pf <C p 2 , 

Eq. (1) can be simplified to p°g = pf + p§ [pi/p^] 1 ^ 2 , so that at low T we have p°g oc T 2 . With increasing T 
this dependence crosses over to p^g- oc exp (— U/2T), i.e. becomes activational. From Eq. (1) we also conclude that 
the effective drag does not follow the evolution of resistivity, p a g, as the classical MIT is continuously swept, due to 
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the variation of the gate voltage. Indeed, the p°g- changes sharply from p" on the metallic side to (piP 2 ) at the 
percolation threshold, and further to p^ on the insulating side. On the other hand, the crossover of p° ff from pf to 
pf is "delayed", as illustrated in Fig. 1. 

The reason why the exact expression for pfg. can be obtained is that the duality transformation [31], which yields 
a closed equation for the p°g- at MIT can be generalized to the case of two layers. This is because, the double-layer 
system can be viewed as a two-component system, in which each component consists of two vertically-separated 
islands, coupled by the mutual drag. Duality transformation between the two types of coupled islands, as depicted in 
Fig. 1, renders a closed matrix equation for the components of the effective resistivity matrix of the two-layer system. 
The corresponding steps arc outlined below. 

Derivation. In the presence of drag, each component of the double-layer system is characterized by its resistivity 
matrix 



pi A I • 2 ~ K pS p'2 

If the two components are equally distributed over the plane, then the effective resistivity matrix, H g, can be found 
exactly. As we demonstrate below the corresponding expression for H e ff has the form 

7? (nnV dI^i + dY 2 ^ 

l-l t .:j 



del 



Dl /2 n! + d\ /2 tz 2 
dI' 2 ^ + D\ /2 n 2 



(p\ /2 + L>2 /2 ) - det (7^ - -fc 2 ) 

where D\ = det(72-i) and D 2 — det(7?-2) are the determinants of the matrices IZi and 1Z 2 , respectively. 
In general, the calculation of the effective resistivity requires the solution of the local Ohm equations 

E a = p a J a + p D J p (4) 
E p = pPJ p + p D J a (5) 

within each double-layer island constituting one of the two components, see Fig. 1. Naturally, Eqs. (4), (5) imply the 
in-plane isotropy of each component. Then it is convenient to view the pairs («7 a , J p ) and (E a , E p ) as two-component 
vectors 



(6) 



and rewrite local equations (4), (5) in the form 



E = 71J, (7) 
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where the matrix, Tl, assumes one of the forms (2) within each component. 

Local equations (7) should be solved together with Maxwell's and continuity equations 



V x E 



_([VxE a ]\_ 
[V x E p ] 



= 



M- 



VJ„ 



= 0. 



(8) 



In order to derive Eq. (3) we demonstrate that, for globally equivalent distributions of the two components, the 
matrix Tl c g- satisfies the following equation 



(9) 



This equation generalizes the Dykhne result [31] to the case of two layers coupled by drag. It is easy to see that in the 
absence of drag, when the matrices Hi and Tl 2 are diagonal, Eq. (9) immediately yields the conventional expressions 
/°eff = an< ^ Pcff ~ V ' P V \P V 2- I n deriving the closed equation (9) for H a \j we follow the line of reasoning put 

forward by Dykhne [31]. Namely, along with J and E, we introduce the auxiliary variables Jd and Ed, defined as 



Jd = A, 



n x E 



E d = A f 



n x J 



(10) 



where Aj and Ae are some constant matrices, and n is the unit vector normal to the layers. It is easy to check that, 
similarly to E and J, the variables, Ed and Jd also satisfy the Maxwell and the continuity equations 



V x E a 







(Vjd) 



0. 



On the other hand, the Ohm's law Eq. (7) dictates the following relation between Ed and Jd 



E„ 



n x J 



n x 



TT^E 



(AeU-'Aj 1 ) J d = TLdJd- 



(11) 



(12) 



At this point we impose the duality conditions. Namely, we require that within the first component Jd and Ed arc 
related via the matrix Tl 2 , and, conversely, within the second component the relation Jd = Tl\Ed holds. If these 
conditions are met, then the equivalent distribution of the first and second components guarantees that on average 
Jd and Ed are related by the same effective resistivity matrix H e g as the average vectors J and E. Quantitatively, 
the duality conditions are expressed as 



Tlx = AeK^Aj 1 , 7^ 2 = AeTL^Aj 1 . 

It is easy to see that these conditions are satisfied by choosing 



Tlx 



Tin 



As a final step, Eq. (9) emerges from the following chain of identities for average fields and currents 



(Ed) 



n x (J) J =A E [nx (n£{E))] = (AeH^Aj 1 ) (J d ) = Tl cS (J d ). 



(13) 



(14) 



(15) 



With Ae — Tli and Aj 1 = TI2, the last identity in Eq. (15) yields Eq. (9). In general, the effective resistivity matrix 
is symmetric, and, thus, is characterized by three unknown elements. As a result, Eq. (9) can be reduced to three 
second-order algebraic equations. It turns our that only two of them are independent. More precisely, the general 
solution of Eq. (9) can be presented in the form 



Jl eS = olTIx +pn 2 , 



(16) 



where a and (3 are the numbers. In order to find these numbers, it is sufficient to derive two relations between them. 
The first relation expresses the fact that the determinants of the l.h.s. and r.h.s. of Eq. (9) are equal. This yields 



dot (a-Ri + pH 2 ) = {DM 112 . 
The second relation emerges upon direct substitution of Eq. (16) into (9) leading to 

oPTlxTl^ 1 + P 2 n 2 Tli 1 = (1 - 2a/3)X, 



(17) 
(18) 
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where X is the unity matrix. It follows from the above relation that nondiagonal elements of the l.h.s. are zero, so 
that 



P 2 



det(K 2 ) _ D 2 



det(72-i) Di 

Upon solving the system Eqs. (17) and (19), we find the following expressions for coefficients a and (3 



a = 



D{ /4 D 3 2 /4 



det |£>2 /2 ^i + D\ /2 n 2 



= 



d\ /4 d\ /4 



det IdI^K! + D\ ,2 n 2 



(19) 



(20) 



Substituting these expressions into Eq. (16), we arrive at the explicit form Eq. (3) of the effective resistivity matrix. 

Applications. In all realistic situations the drag-related nondiagonal components of the matrices (2) are much smaller 
than the diagonal components, which describe the in-plane transport. Under this condition, the effective drag resis- 
tance between the 2D layers can be simplified to 



Peff 



(21) 



The case of drag between a homogeneous layer and a two-component system, considered in the Introduction, corre- 
sponds to p2 = Pi- Then Eq. (21) immediately reduces to Eq. (1). Below we consider two more realizations of the 
double-layer system, in which both layers are strongly inhomogeneous. 

(i) Symmetric layers. This situation (see Fig. 2) emerges when both layers are identical (e.g., positioned symmetrically 
with respect to the donors). Moreover, we will assume for simplicity that the gate voltages applied to the both layers 
are the same. Then, in the vicinity of the classical MIT, the islands (see Fig. 1) will be composed of either two 
metallic or two insulating components. Substituting p\ — p\ and p\ = p% into Eq. (21) we obtain 



1 



PeS- 2 



P2 



1/2 



Pi 



-1/2 



(22) 



In contrast to Eq. (1), pf and p 2 now stand for transresistances between two metals and two insulators. Similar 
to the case of a homogeneous passive layer, outside the MIT region we have p^g- = pf and p^ = p 2 , respectively. 
However, the behavior of p°^ within the transition region is drastically different from that in Fig. 1. Indeed, the first 

term in Eq. (22) contains a small factor [Pi/p^] 1 ^ 2 K cxp(— W/2T), while the second term contains a large factor 
oc exp(W/2T). Thus, despite p% S> pf , at low temperatures the second term will not only dominate, but can exceed 
p 2 . As a result, p°g will exhibit a maximum as a function of V g in the vicinity of MIT, as illustrated in Fig. 2. 
(it) Electron-hole layers. The sign of transresistance in this case is negative [4]. The phenomenon of drag in the 
system of homogeneous electron-hole layers was previously considered in Refs. [5,7,8] with an emphasis on the role 
of interaction-induced correlations between electrons and holes beyond the random-phase approximation. We will 
consider the spatially inhomogeneous situation, assuming that, without disorder, the concentrations of electrons and 
holes are strictly equal. We will also assume that the disorder potential, acting on electrons and holes, is the same. 
The crucial observation is that, due to their opposite charges, electrons and holes "react" differently to the disorder 
potential. The same potential that creates a "metallic lake" of electrons would deplete the corresponding passive 
region of holes, turning them into insulator. As a result, as the MIT is approached, we arrive to the situation, 
depicted in Fig. 3, when the islands consist of pairs of metallic electrons and insulating holes and vice versa. Then, 
substituting p 2 = p° , p\ = p 2 , and p\ — p\ into Eq. (21) we get 



Pl + Pf 



(23) 



It is obvious from Eq. (23) that, since \p°q\ ~ |pf | cxp (-U/2T), the absolute value of the effective drag exhibits a 
minimum near Vg, as illustrated in Fig. 3. 

Discussion. In addition to the drag resistance at MIT, Eq. (3) allows to calculate the drag-induced corrections to the 
effective conductivity of individual layers. For the case of a homogeneous passive layer, considered in the Introduction, 
this correction has a form 



4 



eff v^pf + 

It is noteworthy, that the sign of the correction Eq. (24) is strictly negative. The reason for that is the underlying 
physics of the drag phenomenon. Namely, for electrons in a passive layer, their interaction with electrons in an active 
layer, that gives rise to drag, can be also viewed as an additional source of their scattering. Therefore at the special 
point pi = when the r.h.s. of Eq. (24) turns to zero, it should be expected that the fourth-order correction to 
cr e fj-, neglected in Eq. (24), has also a negative sign. 

Physical explanation of the fact that \p°g\ between the electron-hole layers has a minimum at MIT is straightforward. 
Indeed, when metallic lakes of electrons are located opposite to the insulating regions of holes (see Fig. 3), then, at 
MIT, the current paths in the active layer are perpendicular to those in the passive layer, so that the conditions for 
drag are unfavorable. The origin of maximum of p°^ at MIT for two correlated electron layers, as depicted in Fig. 2, 
is less transparent. One can speculate that the maximum is due to the fact that, at MIT, the current paths in two 
layers are long, and that, due to perfect correlation, each long path in the active layer has its "counterpart" in the 
passive layer. Note finally, that Eq. (3) is exact and takes into account all the orders in p D . Although modeling of 
the classical MIT with two-component mixture is crude, we believe that, due to strong difference in resistances of the 
components, our predictions (1), (22), and (23) for different types of behavior of p°^ across the MIT remain valid for 
realistic situations. 
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FIG. 1. (a) Schematic illustration of the matrix duality transformation. Mutually dual two-layer islands are connected by 
horizontal lines; (b) The resistivity of the active layer (dashed line) and the effective drag (solid line) are depicted as a function 
of the gate voltage for the case when the passive layer is a homogeneous metal. The value p°g at MIT is given by Eq. (1), so 

that p° s /p§ « {pllplf 12 « 1. 
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FIG. 2. The transresistance across the MIT is depicted schematically for two correlated electron layers at low T. The value 
p^g- at MIT is given by Eq. (22). The dependence p^(V 9 ) is the same as in Fig. 1. 
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FIG. 3. The transresistance across the MIT is depicted schematically for electron-hole system. The value p°g- at MIT is 
given by Eq. (23). The dependence Pgff(^s) is the same as in Fig. 1. 
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